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Abstract 

The effect of quadrupole-type surface vibrations on tlie quadrupole response 
function of heavy nuclei is studied by using a model based on the solution of 
the linearized Vlasov equation with moving-surface boundary conditions. By 
using a separable approximation for the residual interaction, an analytical ex- 
pression is obtained for the moving-surface response function. Comparison of 
the fixed- and moving-surface strength functions shows that surface vibrations 
are essential in order to achieve a unified description of the two characteristic 
features of the quadrupole response: the giant resonance and the low-lying 
states. Calculations performed by setting the surface tension equal to zero 
shows that the low-lying strength is strongly affected by the surface tension. 
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I. INTRODUCTION 

There are two systematic features in the quadrupole excitation spectrum of heavy and 
medium- heavy nuclei: the giant quadrupole resonance at an energy huj ~ 63^4" s Me V and 
the lower-energy states that have often been interpreted as surface vibrations. However the 
exact nature of these states is still under debate since they are strongly affected by shell 
effects and this may be taken as an indication that they are not pure surface modes (see e. 
g- [1], P-14). 

There is a vast literature on the quadrupole response of nuclei, here we refer only to the 
nice pedagogical introduction in [2] and to a recent paper [3] where the effect of coupling 
between the motion of individual nucleons and surface oscillations has been studied in a 
model that allows also for non-linear effects and for collisions between nucleons. The scope 
of our present work is more restricted than that of [3] , we aim at studying only the effects of 
coupling between the motion of individual nucleons and surface oscillations of the quadrupole 
type. Our approach is semiclassical and is based on the solution of the linearized Vlasov 
kinetic equation. Solutions of this equation for finite systems have been obtained by using 
different boundary conditions (fixed- and moving- surface [4,5], see also [6]), while the fixed- 
surface solution can give a reasonable picture of the giant quadrupole resonance, it does 
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not satisfactorily describe also the low-lying states, however a unified description can be 
achieved within this model if the moving-surface boundary conditions introduced in [5] are 
employed. Physically, this means that we are including in our model a coupling between 
the motion of individual nucleons and the surface vibrations. In our approach this coupling 
does not involve surface vibrations of different multipolarity. Thus, the present model does 
not include effects that are not taken into account also in the random-phase approximation. 



II. FORMALISM 

A. Fixed-surface solution 

Within the fixed-surface theory, and assuming a simplified residual interaction of sepa- 
rable form, 

v{ri,r2) ^ K,2rfrl, (2.1) 

the quadrupole response function of a spherical nucleus described as a system of A interacting 
nucleons contained within a cavity of radius it! = 1.2^3 fm is given by [7] 

where s = — is a convenient dimensionless variable {vp is the Fermi velocity). The zero- 
order response function 7^-22 (s) is analogous to the single-particle response function of the 
quantum theory and is given explicitly by [8] 

K2is) = 1^- E E ^2iv /' dxx' s^^ix) }^"N{x)r ^2.3) 
Stt ep ^f5f±2 s + ie- s^n (x) 

Here ei? is the Fermi energy, the coefficients are C|o = \ and C|_|_2 = |. The functions 
Snwix) are defined as 

, , nTT -|- arcsin(x) , 

SnN{x) = — (2.4) 

X 

and the quantity £ is a vanishingly small parameter that determines the integration path at 
poles. The Fourier coefficients Q'^{^{x) are the classical limit of the quantal radial matrix 
elements and are given by: 



QIn{x) = + ior{n,N) + (0,0) (2.5) 

^nN\X) ^nN{X) 

Qlix) = R'il - ^x') . (2.6) 

The response function (2.3) involves an infinite sum over n, however in practice it is sufficient 
to include only a few terms around n = in order to fulfill the energy-weighted sum rule with 
good accuracy. The form (2.6) of the coefficient Qoo imphes that the term (n, N) — (0, 0) 
does not contribute to the strength function, hence this term can be omitted from the sum 
in Eq. (2.3). 
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B. Moving-surface solution 



Within the moving-surface theory of [5], the collective response function (2.2) is replaced 

by 

n22{s) ^n22{s)+S22{s), (2.7) 

with 7?.22(s) still given by Eq. (2.2), while S22{s) represents the moving-surface contribution. 
With the simple interaction (2.1) the function 822(8) can be evaluated exphcitly as (see 
Appendix) 

_ [xljs) + K2goR'nl,{s)Y , . 

1 - K2nUs) [C2 - X2{S)][1 - K2nUs)] + ^2R'[X'2(S) + QoR']' ' ^ ' 

with C2 = AaR"^ (a ~ IMcVfm^^ is the surface tension parameter obtained from the mass 
formula) and Qq = A/'^R'^ the equilibrium density. 

The functions X2{^) ^-nd X2{s) are defined as in Refs. [8] and [9] and are given by 

and 

X2{s) — ——€f {s + ie)\- / dxx^ cot[{s + ie)x] -\ — / cixx^f cot[(s -|- i£)a; — 2 arcsin( 
2n ^4: Jo 8 Jo ^ 

+ cot[(s-Fi£)x-F2arcsin(x)])} . (2.10) 

Equation (2.8) is the main result of the present paper, its explicit derivation is lengthy 
but straightforward, the main steps are outlined in the Appendix. Together with Eq. (2.7), 
Eq. (2.8) gives a unified expression for the quadrupole response function, including both the 
high-energy giant resonance and the low-energy excitations. By comparing the two response 
functions (2.2) and (2.7) we can appreciate the effect of the additional surface degree of 
freedom introduced in [5] and in particular the effect of coupling the motion of nucleons 
with surface vibrations of quadrupole type. 



III. RESULTS 



In Fig.l we display the strength function [E — fiuj) 

S{E) ^ --lmTZ{E) , (3.1) 

TT 

obtained for A = 208 using different approximations. The dotted curve is obtained from 
the zero-order response function (2.3) and it is similar to the quantum response evaluated 
in the Hartree-Fock approximation. The dashed curve is obtained from the collective fixed- 
surface response function (2.2). Comparison with the dotted curve clearly shows the effects 
of collectivity: the main strength at about 16 MeV is shifted to lower excitation energy and 
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the giant quadrupolc peak becomes narrower. The strength of the interaction (2.1), chosen 
in order to reproduce the experimental value of the giant quadrupole resonance energy in 
^°^Pb, is «;2 = —1- 10~^MeV fm""^. this value is close to that suggested by the Bohr-Mottelson 
prescription ( [10], p. 509) 

«:2 = -y ^ «^ -0.5 lO-^MeV fm"^ , (3.2) 

(with ujQ = 41^-3 MeV). 

We notice that the width of the giant quadrupole resonance is underestimated by the 
fixed-surface model, this is a well known limit of all mean-field calculations that include 
only Landau damping. Moreover, there is no sign of a low-energy peak in the fixed-surface 
response function. The solid curve instead shows the moving-surface response given by Eqs. 
(2.7)and (2.8). Now a broad bump appears in the low-energy part of the response and a 
narrower peak is situated at the giant resonance energy, thus the moving-surface solution of 
the Vlasov equation introduced in Ref. [5] accounts for both quadrupole modes, although 
only qualitatively. Of course the details of the low-energy excitations are determined by 
quantum effects, nonetheless the present semiclassical approach does reproduce the average 
behaviour of this systematic feature of the quadrupole response. 

Another remarkable feature of the moving-surface response function is that now the giant 
quadrupole peak is narrower than for the fixed-surface solution. This is somewhat surprising 
since we could have expected that introducing a further degree of freedom would result in 
a smearing of the peak, however our result for the giant resonance region is very similar to 
that of the recent random-pase approximation (RPA) calculations of [11] (cf. Fig.5 of [11]). 
Our model does include Landau dampig that, however, turns out to be very small in this 
case, clearly some additional mechanism is required in order to increase the width of the 
giant resonance. Two such mechanisms have been considered in Refs. [12,3], they are the 
coupling to surface vibrations of different multipolarity and the effect of collisions between 
nucleons. It would be interesting to include such effects in the present semiclassical theory, 
however this will be left for future work. 

All the strength functions shown in Fig.l should satisfy the following energy- weighted 
sum rule (EWSR) (see e.g. [10], p. 401): 



dEESiE) = ——AR^ . (3.3) 
47rm 



We have numerically checked that, when integrated up to E' = 30MeV, the response function 
shown by the solid curve exhausts about 98% of this sum rule. The fraction of EWSR 
exhausted by the dashed and dotted strength functions in the same interval is only 80%, 
showing that in these cases there is some more strength at higher energy. 

Another interesting moment of the strength function is the inverse energy- weighted sum 
rule: 



^-.-[.E^-^. (3.4) 

From the three strength functions shown in Fig.l, we have three different inverse moments : 
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n 1 / N 139 1 AR"^ 

m° = — lim7^° s = 3.5 

-1 2 s^o J I40167r ep 

(when evaluating \ims-^oT^22{^) the term {n,N) — (0,0) must be omitted from the sum in 
Eq. (2.3)), 

and 

_ f 140 [g + 2^2771° 

The inverse moment of the zero-order strength function (dotted curve in Fig.l) 
exhausts about 97% of the sum rule (3.5) in the range from to 30 MeV. The collective 
fixed-surface response (dashed curve in Fig.l) instead exhausts almost 99% of the sum rule 
(3.6) in the same energy range while the moving-surface strength function exhusts almost 
100% of the sum rule (3.7), always in the same energy interval. 

It can be seen from Fig.l that, while the fixed-surface response has only one collective 
pole, the moving-surface quadrupole response function displays a two-pole structure. In 
order to get more information about the nature of the low-energy peak, we have performed 
calculations by putting the surface tension parameter a equal to zero. The result is shown 
in Fig. 2, where the dotted curve corresponds to cr = 0. As expected, the giant rcsonacc peak 
is practically unaffected by the surface tension, while the low-energy peak is affected quite 
substantially. The surface tension increases the frequency of the low-enegy peak, which, 
however, is present at a non- vanishing frequency also in the absence of surface tension. In 
the opposite limit, if we let a — > oo, the fixed-surface response is obtained. 

We have performed calculations of the quadrupole response functions also for other values 
of A corresponding to medium-heavy spherical nuclei and the results are qualitatively similar 
to the A = 208 case, so we do not report them here. 

A calculation of the isoscalar quadrupole response similar to the present one has been 
made in Ref. [14], in that case however a rather special external field has been assumed. 
The external force studied in [14] is a pressure that acts only on the surface of the nucleus, 
perhaps this explains why in that case very little strength was found in the region of the 
giant quadrupole resonance. 



IV. SUMMARY AND CONCLUSIONS 



We have obtained an analytical expression for the isoscalar quadrupole response func- 
tion of nuclei that describes qualitatively the two main systematic features of the excitation 
spectrum: the giant resonance and the low-energy states. Our approach is semiclassical and 
is based on the solution of the linearized Vlasov kinetic equation with appropriate boundary 
conditions (moving surface). Comparison of our result (full curve) with the quantum re- 
sponse function of Ref. [11] shows that quantum effects modify substantially the low-energy 
region (a discrete state is obtained in [11] instead of our broad bump), while the giant 
resonance peak is practically the same in the quantum and semiclassical approach. 
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APPENDIX: MOVING-SURFACE RESPONSE 



In our semiclassical approach the nucleus is derscribcd by means of the phase-space 
density /(r, p, t). At equiUbrium the density is /o(r, p) = F(/io(r, p)), where ho is the 
equihbrium mean-field hamiltonian. A weak external driving field T4j;t(r, = (3{t)Q{r) 
induces small fluctuations of the equilibrium distribution /q and of the mean field. Since 
the two fluctuations are related, we have a typical self-consistency problem. For spherical 
system the problem can be solved by means of the partial- wave expansion [4,6] 

5/(r,p,a;)= ^ [5f^\{e,\r,u) + 5f^-^{e,\r,u)]{v^^^^^ (Al) 



LMN 



where ©^^(a, /3, 7) are the Wigner rotation matrices, e = /iq is the particle energy and A 
its angular momentum. 

In the approach of [5] the functions 5fl^^ satisfy the integral equation (see Ref. [6] for 
details) 

g±i[a;r(r)— 7V7(r)] 

SfM^fie^X^r^u) = F\e)— — — -——.ujpr{R)5RLM{^) 



sin[o;r(i?) - N^{R)] 
+{ r dr'Slf^(e, A r')e^y^-<^'^-^^^^'^^ + C^^(e, A, u)}e^^^ 



(A2) 



u!T{r)—Nf(r)] 



with 



B 



L± 
MN 



d 



e,X,r,u;) = F'{e){-± 



iN 



X 



dr Vr{e, X,r) mr'- 



■) [P{u^)QLM{r) + SVLM{r,u;) 



<^ViM(r,.;) = ^^EjW|,|) 
X j dej dXxj ^^-pyVL(r, r')[5f^\f{e,X,r',u;) + 5f^-j^{e,X,r', 



(A3) 
(A4) 



CO 



and 



C^^(6,A,a;) = {e^^i-riR)-N^m T dr5i+v(e. A, r)e-[--W-^^Wl 

which is equivalent to the linearized Vlasov equation with the moving-surface boundary 
conditions 



The surface fluctuations 5Rlm{'^) are related to the functions Sf^j^ by 



5Rlm{<^) 



Stt^ R^ 



LN\ 



2' 2' 



(A5) 



(A6) 



J dej dXXpr{R)[dft-'N{eA,R,uj) + SftrN{e,X,R,uj)-2F'{e)SVLM{R, 



For a separable interaction of the multipole-multipole type, like (2.1), the integral equation 
(A2) can be reduced to an algebraic equation that, in the particular case L = 2 and Vext — 
(3{t)r^Y2M{^), gives Eqs.(2.7) and (2.8) for the quadrupole response function. 
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FIGURES 



FIG. 1. Quadrupole strength function for a hypothetical nucleus of ^ = 208 nucleons. The 
dotted curve shows the zero-order (static moan field) aproximation, the dashed curve instead 
shows the collective response evaluated in the fixed-surface approximation. The full curve gives 
the moving-surface response. 

FIG. 2. The full curve is the same as in Fig.l, the dotted curve has been obtained in the 
moving-surface approach for vanishing surface tension. 
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